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ABSTRACT 

Shipboard  measurements  of  the  vertical  profiles  of  mean 
humidity,  temperature  and  wind  speeds  over  the  ocean  have 
been  investigated  to  verify  their  suitability  to  accepted 
turbulence  similarity  theory,   In  addition,  the  profiles  have 
been  used  to  obtain  estimates  of  the  temperature  structure 
parameter.   One  hundred  eighty-two  profiles  observed  over 
a  seven  month  period  in  1974  aboard  the  R/V  Acania  in 
Monterey  Bay  have  been  compared  to  previous  experimental 
results  obtained  in  defining  the  meteorological  parameters 
of  friction  velocity,  roughness  length,  drag  coefficient 
and  the  Richardson  number.   The  temperature  structure  param- 
eters computed  from  the  profiles  were  provided  to  Lund  (1975) 
for  comparison  with  direct  measurements. 

In  general,  favorable  results  were  obtained  in  the  param- 
eter comparisons  which  provided  additional  support  to  the 
suitability  of  shipboard  profile  measurements  defining  the 
marine  surface  layer.   However,  some  evidence  exists  in 
roughness  lengths  and  drag  coefficient  values  to  indicate 
that  the  ship's  structure  may  be  a  factor  influencing  results 
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I.   INTRODUCTION 

The  air-sea  interface  region  has  many  physical  properties 
yet  to  be  fully  understood.   Only  in  the  past  decade  has 
there  been  a  concerted  effort  to  make  measurements  in  the 
planetary  boundary  layer  over  the  ocean  and  to  examine  exist- 
ing laws  describing  it.   The  Navy  is  actively  engaged  in 
marine  atmospheric  boundary  layer  research  for  many  reasons. 
One  of  which  is  to  accurately  access  the  environmental 
influence  on  optical  wave  propagation  over  the  ocean.   Exper- 
imental evidence  suggests  that  certain  near  surface  atmospheric 
parameters  can  give  reasonable  results  in  specifying  the 
optical  conditions  for  the  first  tens  of  meters  above  the  sea. 
In  this  area,  problems  involving  applied  theory,  shipboard 
instrumentation  and  data  processing  exist. 

A  primary  purpose  of  this  work  is  to  evaluate,  for  the 
marine  boundary  layer,  expressions  which  have  been  formulated 
on  the  basis  of  previous  experiments. 

Data  examined  in  the  study  were  measured  from  the  Naval 
Postgraduate  School  research  ship,  the  R/V  Acania  shown  in 
figure  1.   Beginning  with  accepted  basic  surface  layer  theory 
and  advancing  to  turbulence  similarity  theory,  observed  data 
is  analyzed  in  part  for  its  practicability  in  defining  optical 
wave  properties  in  an  ocean  environment. 

The  suitability  of  the  R/V  Acania  as  a  platform  for  ocean 
micrometeorological  measurements  was  previously  examined  in 
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detail  by  Cavanaugh  (1974)  and  Welsh  (1974)  and  verification 
of  their  work  is  a  secondary  objective. 

Examined  in  this  study  are  ten  minute  time  averaged 
values  of  wind,  temperature  and  relative  humidity  measured 
at  generally  four  levels.   The  observations  were  made  during 
the  months  of  February,  March,  August  and  September  1974 
while  the  R/V  Acania  was  at  anchor  approximately  two  miles 
offshore  in  Monterey  Bay,   Coincident  with  these  measurements 
were  measurements  to  obtain  spectra  of  the  wind  and  temper- 
ature f lunctuations  and  also  measurements  of  optical  propaga- 
tion characteristics.   The  agreement  of  the  profile  data  with 
boundary  layer  theory  predictions  and  results  of  others  will 
be  examined  to  determine  the  validity  of  the  shipboard 
measurements , 

The  boundary  layer  parameters  examined  are  the  roughness 

length  (z  ),  friction  velocity  (u.),  drag  coefficient  (CL.), 
o  *  u 

Richardson  number  (Ri),  and  the  temperature  structure  param- 
eter (C*).   Examinations  of  some  of  these  parameters  are 
accomplished  by  applying  the  data  to  different  calculation 
methods. 

The  structure  parameter  is  calculated  from  an  indirect 
method  based  on  similarity  turbulence  theory  as  proposed 
by  Wyngaard,  Izumi,  and  Collins  (1971).   The  results  were 
used  in  comparisons  by  Lund  (1975),  who  analyzed  direct 
measurements  of  this  parameter  for  the  same  data  periods. 
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II.   THEORETICAL  BACKGROUND 

A .   GENERAL 

Investigations  of  temperature  and  wind  fluctuations  in 
the  surface  layer  provide  a  basis  for  defining  the  stratum 
dynamics  by  turbulence  similarity  theory.   Measurements  of 
mean  wind  (u),  mean  temperature  (T)  and  mean  relative 
humidity  (q)  over  a  time  scale  greater  than  that  of  the 
fluctuating  motion  are  used  to  estimate  the  vertical  fluxes 
of  momentum  and  heat  due  to  the  fluctuations,   This  can  be 
accomplished  because  profiles  of  the  aforementioned  parameters 
are  maintained  by  the  momentum  and  heat  fluxes.   Turbulence 
similarity  theory  which  is  often  referred  to  as  the  Monin- 
Obukhov  theory  relates  the  profiles  to  the  fluctuation 
properties.   Statistical  analysis  is  used  to  arrive  at 
empirical  expressions  which  are  accepted  and  are  consistent. 
Dimensional  analysis  was  used  in  the  formulation  of  all 
expressions. 

The  empirical  relationships  relating  gradients  of  wind, 
temperature,  and  humidity  to  the  two  fluxes  are 

~—   u* 

tfu  _   *  .  ,z.  .1 . 

3i  "  k^  ^ra  <L>  (1) 

3T    T*  .  ,z.  ,0v 

*5  a  e  *h  (l}  (2) 

H  -  E  *e  (!>  (»> 
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Starred  quantities  of  wind  (u+)  and  temperature  (T^) 
include  the  vertical  flux  parameters  as  indicated  in  the 
following  expressions, 

u„  =  (-w)*  (4) 

**  =  -  5  <« 

q*  -  -  5  (6) 


B.   PROFILES  IDENTIFYING  SURFACE  LAYER  CHARACTERISTICS 

Logarithmic  plots  of  mean  wind  profiles  can  be  used  to 
define  certain  surface  layer  properties  such  as  roughness 

length  (z  ),  friction  velocity  (u.)  and  drag  coefficient 

o  * 

(Cn).   The  latter  is  normally  calculated  for  a  ten  meter 
height . 

Integrating  equation  (1)  and  assuming  neutral  or  near 
neutral  stability  conditions,  where  <f>  (z/L)  =  1,  yields  the 
equation 

u*    z 

u(z)  =  —  m  — 

o  v  ' 

The  wind  speed  at  the  surface.  u(z  )  is  assumed  to  be  zero 

'     o 

for  reasons  of  profile  analysis.   However,  that  may  not  be 
the  case  if  the  surface  is  also  in  motion,  as  is  the  case 
over  the  open  ocean,   Wind  measurements  taken  at  several 
levels  and  plotted  as  a  logarithmic  z  profile  yield  an  inter- 
cept of  z   and  a  slope  of  k/u^.   Figure  2  is  a  representation 
of  theoretical  versus  actual  wind  profiles  in  the  marine 
surface  boundary  layer. 


13 


Unstable 


U  Wind    Speed     Im/tec 


Figure   2 


14 


The  drag  coefficient  is  defined  as  the  ratio  of  the  squared 
friction  velocity  to  the  squared  mean  wind  speed  at  a  specific 
height , 

CD  "  -,  <8> 

u 
z 

C.   DESCRIPTIONS  OF  STABILITY  INFLUENCE  ON  THE  SURFACE  LAYER 

The  functions  <t>    ,    d>„,  d)„  are  used  to  account  for  the 

m'   H'   E 

stability  influence  in  equations  (1),  (2),  and  (3),   The 
stability  of  the  surface  layer  regime  plays  an  important 
role  in  determining  flux  intensities,   Monin-Obukhov  theory 
determines  the  level  where  mechanical  and  thermal  production 
of  kinetic  energy  are  equal  to  a  dimensionable  length  (L) , 
The  ratio  of  the  height  of  the  measurement  (z)  to  the  Monin- 
Obukhov  length  (L)  serves  as  a  stability  index,  z/L.   The 
Monin-Obukhov  length  is  defined  as 

T       V    * 

(9) 


where 


T*v  -  T*  +  ,61  q*  T  (10) 


In  unstable  conditions  the  stability  parameter  (z/L)  is 
approximately  equal  to  the  Richardson  number  (Ri)  which  is 
often  used  as  the  stability  parameter,   The  following  relation- 
ships between  the  stability  parameter  and  the  Richardson 
number  have  been  proposed  by  Dyer  and  Hicks  (1970)  (unstable) 
and  Webb  (1970)  (stable). 
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z 

] 

Rl 

L 

1    - 

a   Ri 

z 

L 

Ri 

(stable)         (11) 

(unstable)       (12) 

a  is  an  empirically  derived  constant  equal  to  ,5, 

D.   INDIRECT  CALCULATION  OF  THE  TEMPERATURE  STRUCTURE 
PARAMETER 

The  temperature  structure  parameter  (C*)  is  the  most 

important  variable  in  determining  the  optical  propagation 

characteristics  of  the  atmosphere,   One  means  for  calculating 

C£,  is  by  employing  turbulence  similarity  theory.   Since 

atmospheric  refractive  index  fluctuations  are  not  included 

in  these  calculations,  the  method  indicated  is  an  indirect 

estimation  of  atmospheric  optical  properties.   In  similarity 

theory  estimates  of  the  dissipation  of  turbulent  kinetic 

energy  (e)  and  dissipation  of  temperature  variance  (efi) 

provide  the  link  between  direct  methods  such  as  proposed  by 

Corrsin  (1951),  equation  13,  and  the  conveniently  relatable 

vertical  fluctuations  of  wind  and  temperature  in  the  surface 

layer, 

<f>T(ki)  =  B^eeg  ~1/3  k~5/3  (13) 

Research  by  Wyngaard,  Izumi  and  Collins  (1971)  has  resulted 
in  nondimensionalized  expressions  relating  atmospheric 
stability  to  the  structure  parameter.   Dimensional  analysis 
of  the  height  (z),  the  friction  velocity  (u+),  the  mean 
potential  temperature  gradient  (-*—),  e,  en  and  Cz,   produce 

oZ  u         1 
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functional  expressions  dependent  only  on  the  Richardson 
number.   Similarity  principals  are  elucidated  by  using  the 
nondimensional  stability  parameter  Richardson  number  defined 
as  equation  14  for  indirectly  obtaining  a  comparable  structure 
parameter  to  identify  the  optic  properties  existing  in  the 
surface  layer. 


Ri  = 


(lu}2 


(14) 


Equations  including  the  Richardson  number  functions 
(fi,  f2,  f3)  are 


£|  =  f   (Ri)  (15) 


=  f,  (Ri)  (16) 


-*(H)2 


'T    _  3.2  f2 


C2 
J7~  =  *'*k   l2  =  fa  (Ri)  (17) 

zh/3  (|6  2    f  J  2 

The  latter  expression  relates  the  temperature  structure 
parameter  and  stability  parameter  and  allows  for  an  empirical 
determination  of  the  f3  function  at  varying  conditions  of 
stability.   If  the  function  f3  is  defined,  the  temperature 
structure  parameter  can  be  readily  calculated  from  more  easily 


measured  parameters,  -^—  and  Ri 

r  dZ 
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III.   EXPERIMENTAL  OBSERVATIONS 

A.  DATA  COLLECTION 

The  Naval  Postgraduate  School  R/V  Acania  was  the  obser- 
vational platform,   The  ship  was  anchored  in  Monterey  Bay 
at  a  location  approximated  in  figure  3,   This  location 
facilitated  complimentary  and  simultaneous  optical  experi- 
ments from  two  shore  stations,  Hopkins  Marine  Station  and 
the  Holiday  Inn.   Normally,  observational  periods  were  three 
days  in  duration  beginning  with  a  one  day  equipment  calibra- 
tion and  set  up,  continuing  with  a  one  day  on  station  systems 
test  and  finishing  with  one  day  of  data  collection,   Data 
collection  was  mainly  during  evening  hours.   The  entire 
operation  required  considerable  inter-departmental  coordination 
between  the  meteorology  and  physics  departments. 

Data  was  collected  and  analyzed  for  the  dates  and  time 
shown  in  table  1.   These  observations  represent  over  ninety 
percent  of  all  the  mean  profiles  obtained  in  Monterey  Bay 
during  the  period  from  25  February  to  20  September  1974, 

Figure  4  shows  the  general  instrument  configuration 
used  during  the  experiments,   Variations  in  the  arrangement 
occurred  primarily  because  of  individual  component  failures, 

B,  EQUIPMENT 

The  mean  wind  measurements  were  made  with  a  Thornthwaite 
Associates  cup  anemometer  wind  profile  register  system, 
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DATE 

TIME 

NUMBER 
OF 

NUMBER 

OF 

LEVELS 

PROFILES 

26  February 

1113-2225 

4 

23 

27  February 

1643-2150 

4 

25 

14  March 

1703-2138 

2 

20 

27  March 

1823-2142 

3 

17 

28  March 

1824-2150 

4 

14 

13  August 

1630-2200 

4 

26 

18  September 

1617-2202 

4 

24 

19  September 

1724-2353 

4 

25 

20  September 

0004-0644 

4 

8 

Table    I 
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Level  Variable 


G  Tq 


GTq 
uf  q 

u  T  q 


To 


-^ 


Sea    Surface 


Height  (meters) 


--13.9 


7.6 
6.6  0 


■-4.2 


-0 


Figure   4 
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model  number  104.   In  operation  the  shaft  of  a  three  cup 
anemometer  unit  serves  as  the  shutter  between  a  light  source 
and  a  photocell  for  each  revolution.  The  cups  are  plastic 
cones  reinforced  with  aluminum  frames.   They  are  attached 
to  the  rotating  shaft  by  stainless  steel  tubing  spaced  at 
120  degree  intervals  about  the  shaft  as  shown  in  figure  5. 
The  three  cup  assembly  sets  along  with  the  other  sensors 
were  positioned  at  four  levels  on  the  bow  with  electrical 
leads  to  the  after  deck  house  laboratory.   The  sets  have  the 
characteristics  of  low  starting . speeds  with  a  small  amount 
of  internal  friction  which  aids  in  checking  inertial  over- 
shoot. 

Temperature  sensitive  quartz  crystal  probes,  (Hewlett 
Packard  model  HP-2850)  were  used  to  measure  mean  temperatures 
at  each  level,   RF  signals  from  the  crystal  probes  and  from 
a  reference  oscillator  were  mixed  in  the  HP-2801A  readout 
unit  to  produce  a  beat  frequency  whose  signature  can  be 
analyzed  to  within  0.001  degrees  centigrade  per  hertz. 
Each  sensor  simultaneously  received  pre-experiment  calibra- 
tion against  a  platinum  resistence  wire  thermometer  in  a 
temperature  controlled  circulating  water  bath  over  the 
expected  temperature  range.   The  accuracy  in  achieving  a 
.005  degree  centigrade  correction  factor  was  a  constant  for 
each  probe.   A  3.7  meter  flexible  coaxial  cable  is  permanently 
attached  to  the  sensor  head  and  the  mast  mounted  probes  are 
housed  in  an  aspirated  shelter  as  depicted  in  figure  6, 
Temperature  values  were  automatically  recorded  on  a  printer 
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Figure  5 
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C.C.  Breidert  Company 
Air-X-Houster  Type  6L 


Dunmore-type  Lithium 
Chloride  Sensor 


Quartz  Thermometer  Probe 


Figure  6 
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tape,   The  Hygrodynamics  Digital  I  hygrometer  indicator  and 
the  Dunmore-type  lithium  chloride  sensor  were  used  to  measure 
the  relative  humidity,   The  Dunmore-type  lithium  chloride 
sensor  is  illustrated  in  figure  6,   The  equipment  operates 
on  the  basis  of  resistence  change  within  an  electrolytic 
solution  generating  a  reference  voltage  variance  which  is 
proportional  to  the  relative  humidity  change.   Automatic 
temperature  compensation  in  the  instruments  meet  the  follow- 
ing specifications  for  relative  humidity, 

+   3%  relative  humidity  below  90% 
+   4%  relative  humidity  above  90% 
Calibration  of  the  sensors  was  accomplished  by  comparative 
methods  using  saturated  salt  solutions  in  an  enclosed 
chamber. 

Relative  humidity  was  observed  and  recorded  similar  to 
that  of  temperature,  as  a  printed  output.   A  more  detailed 
description  of  the  devices,  their  calibration  and  deploy- 
ment can  be  obtained  from  Cavanaugh  (1974), 
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IV,   ANALYTICAL  DATA  PROCESSING 

A.   DATA  EDITING 

During  the  experiment,  mean  wind  measurements  over  ten 
minute  intervals  were  recorded  in  a  log  book,   The  recordings 
were  the  number  of  cup  revolutions  for  the  interval  as 
recorded  by  electronic  counters.   Temperature  readings  were 
printed  on  a  paper  tape  from  a  monitoring  digital  readout 
at  approximately  one  second  intervals,   Initially  relative 
humidity  measurements  were  recorded  the  same  as  the  tem- 
peratures.  However,  due  to  a  printer  malfunction,  these 
readings  were  recorded  by  hand  in  log  book  after  the  28 
March  cruise.   All  information  was  later  coded  and  key 
punched  on  computer  cards  for  further  analyses, 

Preliminary  screening  processes  were  performed  on 
computer  listings,   The  data  were  examined  for  gross  errors 
or  inconsistencies,   Very  few  deletions  were  necessary  at 
this  phase.   Next,  individual  mean  wind  and  mean  temperature 
measurements  were  plotted  to  give  logarithmic  graphical 
profiles  for  each  time  period.   The  computer  program  devised 
for  this  operation  also  performed  a   least  squares  fit  of  the 
data. 

In  addition,  mean  wind  profiles  for  certain  periods 
were  closely  examined  for  any  biases  which  might  be  caused 
by  the  ships'  hull  or  superstructure.   Approximately  ten  per 
cent  of  all  the  measurements  was  discarded  either  because 
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of  incongruous  profiles  or  uncomplementary  data  times. 
Rejected  profiles  were  mainly  those  of  wind.   They  occurred 
principally  when  winds  were  very  light. 

Temperature  profiles  were  rejected  when  one  or  more  of 
the  automatically  recorded  readings  appeared  to  be  spurious. 
This  was  not  discovered  during  the  initial  data  inspection. 

B.  DATA  RECORDING  ADJUSTMENTS 

Since  the  data  was  obtained  over  a  seven  month  period, 
significant  alterations  in  data  recording  took  place  due  to 
equipment  failure  or  experiment  deviations.   Mean  wind  read- 
ings were  obtained  at  four  levels  for  all  the  periods  except 
when  the  number  of  levels  were  two  and  three  for  March  14 
and  27  respectively.   In  August,  the  order  in  which  the  mean 
wind  levels  were  labeled  changed  from  previous  cruises.   This 
data  variation  was  compensated  for  by  analyzing  the  measure- 
ments on  a  cruise  to  cruise  basis  with  simple  computer 
program  modifications.   There  were  no  major  problems  encoun- 
tered in  programmed  analysis  of  the  temperature  data,   When 
the  indicated  relative  humidity  was  missing,  adjustments  were 
made  by  taking  an  average  between  the  two  adjacent  available 
measurements . 

C.  PROGRAMMING  CONSIDERATIONS 

Cup  anemometer  readings  were  recorded  in  revolutions  per 
time  period,   The  relationship  between  revolutions  to  speeds 
in  nonlinear.   It  changes  with  wind  speed,   A  computer 
subroutine  was  devised  to  convert  the  number  of  cup  rotations 
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into  wind  speed.   The  conversion  accuracy  for  wind  speeds 
decreases  with  decreasing  revolutions, 

An  IBM  System/360  source  library  program  called  LSQPL2 , 
Least  Squares  Polynomial  Fitting,  was  an  integral  part  of  the 
data  analyses.   LSQPS2  not  only  served  to  give  linear  func- 
tions for  data  profiles,  but  also  was  used  as  a  means  for 
approximating  roughness  length  (z  )  values,   Once  the  least 
squares  linear  profiles  had  been  calculated,  profile  values 
were  estimated  for  levels  one  meter  either  side  of  the 
observation  levels.   This  was  done  for  the  purpose  of  defin- 
ing gradients  used  in  computations  of  the  Richardson  numbers. 

The  basis  for  describing  stability  in  discerning  optical 
propagation  conditions  is  the  Richardson  number.   In  testing 
for  the  feasibility  of  applying  identical  relationships  over 
the  ocean,  Richardson  numbers  were  determined  from  the 
Monterey  data  utilizing  procedures  similar  to  those  set 
forth  by  Wyngaard  et  al,  (1971),   Mean  profiles  of  temperature 
and  wind  were  obtained  by  an  analytical  least  squares  fit  of 
ten  minute  mean  observations  at  their  respective  instrument 
heights.   This  provided  the  wind  and  potential  temperature 
gradients  necessary  to  compute  a  Richardson  number  at  each  of 
the  four  instrument  levels.   Equation  18  is  the  approximation 
for  the  mean  wind  speed  gradient  between  two  levels  used  in 
equation  14, 
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3u 

dz 


=  &}*  =  u(z2 )  ~  u(zt  ) 

Az    z3  in  z2/Zj 


(18) 


z3  represents  the  geometric  mean  height,  (z1    •  z2)2,   A 
similar  expression  was  used  to  compute  temperature  gradient, 
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V.   DATA  ANALYZATION  AND  COMPARISONS 

A.   FRICTION  VELOCITY  AND  ROUGHNESS  LENGTH  RESULTS 

Important  parameters  describing  the  marine  boundary 
layer  are  the  friction  velocity  (u^)  and  roughness  length 
(z  ).   The  friction  velocity  is  a  measure  of  the  momentum 
transfer  described  by  equation  4  and  the  roughness  length  is 
a  boundary  parameter  defined  as  the  height  where  the  loga- 
rithmic profile  becomes  a  valid  specification.   These  param- 
eters were  calculated  from  the  data  by  two  methods.   One 
method  utilized  all  profile  levels,  the  other  utilized  only 
the  wind  speed  at  one  level  and  the  air  water  temperature 
difference.   This  was  done  to  determine  what  difference  would 
occur  if  the  surface  layer  properties  were  obtained  using 
data  from  one  level  instead  of  from  a  four  level  profile. 

The  first  method  was  based  on  a  least  squares  fit  of  a 
stability  corrected  four  level  profile.   The  stability 
influence  was  determined  at  each  level  from  the  Richardson 
number.   The  following  Richardson  number  categories  were 
used  in  the  corrections. 

-0,05  <  Ri  <   ,001        neutral 
Ri  >   .001        stable 
Ri  <  -,05         unstable 
Stability  relationships  in  equation  11  and  12  were  used  to 
obtain  the  stability  length  (z/L)  which  was  then  used  to 
obtain  the  correction  to  the  empirical  expression  for  mean  wind. 

30 


The  plot  of  the  mean  wind  versus  the  logarithm  of 
height  with  the  stability  correction  should  theoretically 
give  a  straight  line  with  a  slope  of  k/u  +  and  an  intercept 
of  z  ,   Equation  19  was  used  to  incorporate  the  stability 
influence  into  the  profile 

in  z  -  ♦  C2-)  =  -  u(z)  +  In  z  (19) 

x       YmvL/    u^   v  '  o  v  ' 

The  function  iKy-)  was  defined  as  follows  from  Paulson  (1970) 

♦  (f)  =  2  in  p^P]+  inp^]  -  2tan~1  x  +  \      (20) 

x  =  (1  -  16(f))*        unstable  (21) 

♦(f)  =  -4,7  f  stable  (22) 

A  least  squares  fit  was  applied  to  the  corrected  profile  to 
obtain  the  best  estimate  of  the  slope  and  intercept  thus 
yielding  the  profile  values  for  friction  velocity  and  rough- 
ness length, 

The  second  method  was  to  use  a  calculation  suggested  by 
Cardone  (1969).   He  formulated  an  iterative  procedure  to 
find  friction  velocity  and  roughness  length  given  a  mean 
wind,  an  air-sea  temperature  difference  and  the  instrument 
heights.   By  utilizing  wind  and  temperature  gradient  expres- 
sions combined  with  equations  similar  to  9,  an  iterative 
algorithm  was  developed,  which  gives  stability  length 
convergence  to  a  desired  accuracy.   The  friction  velocity 
was  obtained  using  the  stability  criterion.   The  roughness 
length  was  then  calculated  from  the  friction  velocity  and 
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the  following  expression  from  Cardone, 

zo  =  n;  +  cau*  -  c3  (23) 

The  constants  contained  in  equation  23  were  derived  so  that 
the  ten  meter  drag  coefficient  (C1f))  would  be  a  minimum 
when  the  mean  wind  speed  was  approximately  six  meters  per 
second  under  neutral  conditions. 

The  results  from  the  two  methods  are  compared  in  figure 
7  which  is  a  logarithmic  plot  of  z   versus  u+.  Comparison 
is  also  made  with  Charnock's  (1955)  expression. 


u*2 


z   = 


o    a  g  a  =  81,1  (24) 

The  similar  curvature  of  all  three  plots  indicates  agree- 
ment on  the  shape  of  the  friction  velocity  and  roughness 
length  relationship,   However,  the  least  squares  method 
produced  roughness  lengths  on  the  order  of  one  magnitude 
larger  than  Cardone ' s  values, 

In  general,  both  methods  appear  to  give  good  results, 
however,  the  average  values  obtained  by  the  least  squares 
method  show  less  dispersion  at  low  friction  velocities  than 
does  Cardone ' s  one  level  method, 

A  comparison  of  the  roughness  length  versus  the  ten 
meter  mean  wind  speed  was  made  between  results  obtained  by 
Ruggles  (1970)  and  the  least  squares  method  as  shown  in 
figure  8.   Ruggles  hypothesised  that  the  z   peaks  were 
associated  with  the  effects  of  wind  wave  coupling.   He 
designated  the  first  peak  as  being  coincident  with  the 
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minimum  wind  necessary  for  wave  generation,  The  results  of 
this  study  show  similar  maximums  and  minimums,  most  notibly 
the  maximum  near  3.5  meters  and  minimum  near  4.5  meters. 

B.   DRAG  COEFFICIENT  RESULTS 

Ten  meter  drag  coefficients  (C-iq)  were  calculated  to 
provide  further  verification  of  the  compatibility  of  the 
experimental  profile  results  to  results  of  others,   Friction 
velocities  which  were  computed  by  the  least  squares  method 
and  also  using  Cardone's  expression  and  the  ten  meter  mean 
wind  (uir))  were  used  to  obtain  C.  n  by  the  expression 

cio  =  502  (25) 

Plots  of  C.  n  versus  u^  and  C.  n  versus  Richardson  number 
appear  in  figures  9  and  10.  Superimposed  on  the  graphs  are 
the  results  of  others. 

In  general,  some  aspects  of  the  results  compare  favor- 
ably with  results  of  others  and  indicate  that  over  ocean 
shipboard  application  of  turbulence  similarity  theory  might 
be  valid  and  feasible  in  the  marine  surface  boundary  layer. 
The  least  squares  method  produced  average  drag  coefficients 
that  were  higher  than  the  values  of  other  experiments.   This 
is  attributed  to  the  high  values  of  friction  velocity  possibly 
due  to  the  ship's  influence,   The  problem  was  addressed  by 
Cavanaugh  (1974),  however,  his  results  were  not  as  pronounced 
as  those  in  this  study, 
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C,   TEMPERATURE  STRUCTURE  PARAMETER 

The  most  important  element  for  parameterizing  optical 
propagation  conditions  in  the  marine  surface  boundary  layer 
is  the  temperature  structure  parameter  (C2),   The  indirect 
method  was  adopted  to  define  C2  at  the  four  instrument  levels 
where  a  computed  Richardson  number  was  available,   First,  it 
was  necessary  to  use  the  Richardson  numbers  to  define  the 
f3  function  in  equation  17  that  had  been  empirically  derived 
from  the  Kansas  field  data.   Equation  17  was  then  used  to 
obtain  C2,   Subsequently,  the  temperature  structure  parameter 
was  transformed  into  the  refractive  index  structure  parameter, 
C2  for  additional  comparisons  with  theoretical  predictions, 
spectral  and  optical  results, 

Calculated  C2  decrease  with  height  as  expected  for  all 
stability  conditions.   Above  a  Richardson  number  of  ,20,  the 
f3  function  in  equation  17  is  undefined, 

Temperature  structure  parameter  (C,^)  results  were 
compared  to  direct  measurement  values  by  Lund  (1975), 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  ship's  motion  and  structure  appeared  to  have  had  an 
effect  on  the  wind  profiles.   This  conclusion  is  based  on  the 

fact  that  the  values  of  z   and  C.  ~  obtained  by  a  least  squares 

o      10  J  * 

fit  profile  data  were  consistently  higher  than  values  observed 

in  other  experiments. 

Cavanaugh  (1974)  observed  similar  results  and  alluded 

to  ship's  motion  influence  on  a  logarithmic  profile.   Hull 

interference  could  disrupt  the  wind  stream  passing  the 

instruments.   This  obstruction  could  cause  the  lower  wind 

cups  to  slow  and  thus  increase  the  apparent  gradient.   This 

would  result  in  higher  z   values.   This  possible  influence 

b     o  * 

was  examined  in  the  analysis  by  disregarding  the  lower  sensor 
level  in  the  regression  calculations.   Roughness  lengths  were 
lowered  slightly  in  most  profiles  examined,  however,  not 
enough  to  justify  dropping  this  level  in  profile  analyses. 

A  suggestion  for  future  experiments  would  be  to  position 
a  sonic  anemometer  vertically  over  the  side  of  the  hull  at 
main  deck  level  to  measure  the  forced  vertical  air  motion 
caused  by  hull  wind  blockage.   From  data  obtained  in  this 
manner,  a  ships  wind  stream  signature  for  varying  wind  speeds 
might  be  found.   This  could  find  application  in  the  more 
complex  study  of  measuring  turbulent  motion  while  underway. 

With  regard  to  temperature  structure  parameter  estimates, 
the  function  f   is  required.   Sufficient  data  is  now  becoming 
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available  to  permit  evaluation  of  equation  17  for  over  the 
ocean  observations.   Comparison  of  the  Ci,  spectral  results 
obtained  by  Lund  (1975)  and  those  computed  in  this  study 
could  accomplish  part  of  this  task. 

Because  of  the  quantity  of  data  to  be  screened  only  a 
cursory  analysis  of  the  observations  has  been  completed  in 
this  study.   More  detailed  examinations  into  the  measurements 
would  certainly  yield  additional  knowledge,  particularly  with 
regard  to  the  ship's  motion  and  structure  influences  and 
instrument  biases. 

It  would  be  desirable  to  have  mean  wind  observations 
closer  to  the  ocean  surface.   Combined  with  the  present 
profiles,  a  more  detailed  description  of  wind-wave  coupling 
and  its  effect  on  turbulent  flow  could  be  undertaken,   The 
additional  measurements  could  be  made  from  a  buoy  positioned 
upwind  of  the  ship.   Another-  recommendation  which  is  unrelated 
to  application  of  turbulence  similarity  theory   is  that  a  small 
portable  computer  be  obtained  for  future  experiments  to 
facilitate  real  time  data  analysis.   The  system  now  employed 
is  extremely  inefficient  because  many  data  discrepancies  are 
not  recognized  until  after  the  cruise, 
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